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I. INTRODUCTION 



Systems that are mean-field or near-mean-field are common in nature. Examples of such 



systems include metals with long-range elastic forces 



rr 



, |2| , earthquake faults with long-range 
stress transfer Green's functions [3|, and polymers [4?]. The connection between the range of 
the interaction and mean-field behavior was made by Kac and collaborators [5|] who noted 
that a system with a pairwise additive potential of the form 



V{x) = Vr{x)+^''^{^x), 



(1) 



becomes mean-field in the limit 7^0. In Eq. ([T]) x = |x|, Vr(x) is a short range reference 
potential, and d is the spatial dimension. The limit 7 — is taken after the thermodynamic 
limit and before a critical point is approached. It is also required that 



j rfx7'^|$(7x)| < 00. 



(2) 



so that the energy per particle or spin remains finite in the 7^0 limit. The interaction 
range R is defined by the second moment of the potential, 



(3) 



Hence, as 7 — > 0, i? ^ 00. We will refer to systems with R ^ 1 but not infinite as 
near-mean- field] systems with i? — > 00 are mean-field jsl. 

The kinetics of phase transitions is different in systems with R^ 1 than in systems with 
i? ~ 1. For example, nucleation in the former often occurs near a pseudospinodal 0, 0, ^ 



where the sur: 
structure 

R ~ : 



ace tension is small, which results in a nucleating droplet that has a different 



11 



19|. 



12, 



13 



14 



15 



16 



171] than that near the coexistence curve in systems with 



In addition, the early stage growth of the peak of the equal time structure function during 
continuous ordering and spinodal decomposition in systems with R^ 1 is described by the 



Cahn-Hilhard-Cook (CHC) theory 



20 



21 



22l | for a time proportional to InR after the 



quench 14]. The morphology of the early stage evolution differs from that in systems with 
i? ~ 1 j23l . 2^ . 25 ] for which there is no time interval when the CHC theory is applicable jiol ]. 



* Formerly at Department of Physics, Clark University, Worcester, Massachusetts 01610, USA 
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The mean-field limit of several earthquake fault models can be described by an equilibrium 



theory 



26 



271]. In near- mean- field systems the smaller earthquake events are related to 



fluctuations about the free energy minimum near the pseudospinodal 3|. 

In these and other examples the structure of the fluctuations near the mean-field critical 
point and the pseudospinodal is important for understanding the behavior of the system. 
In this paper we analyze the structure of the fluctuations and its relation to the underlying 
clusters. We use field theory, scaling arguments, and cluster analysis and relate the structure 
of the fluctuations to the nature of nucleation, the possible existence of a pseudospinodal 
in supercooled fluids, and the behavior of the models of earthquake faults. The results of 
simulations done to test the predictions are also discussed. 

In Sees. HI] and IIIII we discuss the Landau- Ginzburg theory {3] and the Parisi- 



Sourlas [29|, |30|] approach based on the Langevin equation with random Gaussian noise to 
study fluctuations near mean-field critical points. We use the same field theory techniques 
in Sec. HV] to discuss fluctuations near the spinodal. In Sec. |V] we discuss the fluctuation 
morphology for mean-field and near-mean-field systems. In Sec. I VI I we use the Landau- 
Ginzburg and Par isi- Sour las approaches to discuss the relation of the fluctuations to the 
clusters. We examine the relation between the fluctuation structure and spinodals in su- 
percooled fluids in Sec. IVIII and discuss the relation between the fluctuation structure and 
nucleation in Sec. IVIIII In Sec. IIXI we relate the cluster structure to cellular automata 
models of earthquakes. We summarize our results and discuss future work in Sec. |Xl The 
mapping of thermal systems onto percolation models is discussed in Appendix |Al 

Our main results include the following. (1) There exist objects, which we call funda- 
mental clusters, that have a density and lifetime dependence that is very different from 
the scaling of the density and lifetime of the fluctuations. The fundamental clusters are 
defined by the mapping of the critical point (or spinodal) onto a percolation transition. 
This difference is in contrast to non-mean-field systems where the clusters are geometrical 
realizations of the fluctuations. (2) The fundamental clusters are physical objects that have 
measurable consequences, which are explored for earthquake fault models, nucleation, and 
the measurement of the (pseudo) spinodal in near-mean-fleld systems. 
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II. SCALING OF ORDER PARAMETER FLUCTUATIONS 



We first discuss (near-)mean-field systems 



the Landau- Ginzburg- Wilson Hamiltonian 



Tom the perspective of field theory based on 



31 



32| 



[V0(x)]' + e0'(x) + 0^(x) - /i0(x) 



(4) 



Without loss of generality we have set the proportionality constant in Eq. ([3]) equal to one. 
The partition function Z is 



6(f) e 



^I3H{4>{^)) 



and the probability of the order parameter density 



X IS 



(5) 



(6) 



where j3 = {ksT)'^, T is the absolute temperature, and ks is Boltzmann's constant. 

Because we are interested in (near-) mean-field systems, we scale all lengths with R. We 
first discuss the critical point and defer our discussion of the spinodal to Sec. [IVl To study 
the critical point we set /i = and assume that 

T-T, 



< 1, 



(7) 



where is the critical temperature. For e <^ 1 and h = we can use scaling arguments. It 
is straightforward to see from Eq. (jlj) that 



HU) = rW'-''/' I dy [V0(y)] ' ± 4>\y) + 4>\y) 



(8) 



where y = x/i^e"^/^, 0(x) = e~^^'^(f){x), V = -RV, and the -|- (— ) sign corresponds to e > 
(e < 0). We take e > in this section except where otherwise noted. 

The integral in Eq. ([5]) can be evaluated using saddle point techniques for /^<^g2-d/2 ^ 
We can give this requirement a physical meaning from the Ginzburg criterion [28!], which 
states that a system can be considered to be mean-field if the mean square fluctuations of 
the order parameter are small compared to the square of the order parameter 321]. The 
order parameter is given by L'"^] c/x0(x), where L is the linear dimension of the system 
and (p corresponds to the magnetization in the Ising model. 

The correlation length ^ is proportional to the linear spatial extent of the order parameter 
fluctuations. The mean square fluctuations in the order parameter are characterized by ^"'x. 
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where x is the isothermal susceptibihty 32|. Near a mean- field critical point we have 33 1 

e ~ Re-^'^ (9a) 
0~|e|i/2 (e<0) (9b) 
X-t-\ (9c) 

(Equation ( I9b[) is derived following Eq. (l56l).) The Ginzburg criterion requires that 

0. (10) 

If we substitute the scaling forms in Eq. ^ into Eq. flTUl) . we obtain 

G = ^ (11) 

We will refer to G = /^'^g2-(i/2 ^j^g Ginzburg parameter. In the limit G ^ oo the system 
is mean-field. The system is near- mean- field for G ^ 1 (but finite). The latter criterion 
implies the well known result that the upper critical dimension at the critical point above 



which the system has mean-field critical exponents for all i?, including i? ~ 1, is four 
From Eqs. ([6]) and ([8]) we have 

o _ exp {- m2-'^/Vrfy[i[V0(y)]2 + ^^(y) + 04(y)]| 



32|. 



(12) 



For G = R^e^ "'/^ ^ 1 the Hamiltonian in Eq. (jl]) can be approximated by a Gaussian 32| : 

^g(0) = ^ , (13) 

where Zq is the functional integral over of the numerator in Eq. (fT3ll . We use Pg{4') to 
calculate the structure function S{k): 

e / 50(k) exp[-/5i?^e2-^/V ^^(^2 ^ i)0(k)0(-k)]0(k)0(-k) 



^(/t) = e<0(k)0(-k)> 



(14) 



where k = i?e ^/^k and 0(k) = 0(-k). For e > and h = 0, <0(x)> = 0. We have 



The Fourier transform of Eq. f[T5|) gives the pair distribution function, which we write in 
terms of unsealed variables: 

o^^Hx) - - e"''/^ (16) 
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The {x/^Y~'^ dependence in Eq. (fT6l) is valid for c? > 3. For d = 2 this dependence is 
replaced by in ci = 1 there is no x-dependence in the denominator. The integral 

fdxp^'^\x) is proportional to [e/ R'^e^^'^^'^)^'^ for all d. For scaling purposes we can treat 
p*^^''(x) for X ^ ^ as a constant. Because p*^^^(x ^ ^) is proportional to the square of the 
density of a fluctuation, we see that the fluctuations in the order parameter density scale as 

0f(x ^ i) ^ ^j^d^2-d/2y/2 ^ GiJ^' ^^'^'^ 

Note that the density of a critical phenomena fluctuation does not scale as e^^'^ as might be 
expected from a simple extension of how the order parameter scales for e < in Eq. (]9bp . 
We will discuss this point more fully in Sec. IIVI The scaling of 0f (x) with G^^^^ in Eq. f|T7|) 
justifies the neglect of the 0"^ term in Eq. (fT2!) and the Gaussian approximation in Eq. (fT3ll . 
The susceptibility x is related to the pair distribution function p^'^'> by 32, 



Xocy"dx[p(2)(x)-02], (18) 

If we use the scaling form f lTTl) of 0f(x) in Eq. ffTSjl . that is, x ~ we find x ~ 

consistent with Eq. (!9c|) . 

To show that 0f and have similar scaling behavior in a system with i? ~ 1 and d < 4, 
we again assume that p^^^ (x ^ Qis a constant so that we can write Eq. (|T8l) as x ~ 0/^^^ and 



~ e '^e^'^. Hyperscaling 



32 



33| gives 7 + 2/3 = dz/ so that (p'j ~ e^^ and hence 0j ~ €■ 



/3 



We next discuss the Landau-Ginzburg and Cahn-Hilliard-Cook equations in (near-)mean- 
field systems. We can obtain these equations by noting that the time rate of change of 
an order parameter such as the density is related to the chemical potential p. If the or- 
der parameter is not conserved and there are no other conservation laws (model A in the 
Hohenberg-Halperin classification scheme [m]), then 

oc -p and p oc , (19) 

where F{(f)) is the free energy. We take F{(f)) to be equal to the Landau-Ginzburg- Wilson 
Hamiltonian in Eq. (jl]), which is correct for mean-field systems and a good approximation 
in near-mean-field systems, and assume that the relations in Eq. ( |T9l) are valid in a spatial 
and time dependent context and that the functional derivatives are with respect to 0(x, t). 



19 



28|: 



In this way we obtain the Landau-Ginzburg equation 

= -M^[-i?2v20(x, t) + 2e0(x, t) + 403(x, t) - h] + r]{^, t), (20) 
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where we have added a noise term ?7(x, t). In the remainder of this section e can be either 

n 

positive or negative. For a conserved order parameter (model B ^]) we have 



dt 



oc V ■ J and J oc V/i(x, t). (21) 



If we again interpret the right-hand side of Eq 



we obtain the Cahn-Hilliard-Cook equation 19(] 



jH) as a free energy and include a noise term, 



= MBV=^[-i?W(x, t) + 2e0(x, t) + 40='(x, t) - /i] + r/,(x, t). (22) 

The quantities Ma and in Eqs. (120|) and fl22l) are mobilities and will be discussed in 
Sec. [ml 

To obtain Eqs. (|T9l) and ( 12T1) we assumed local equilibrium; that is, within the coarse 
grained volume used to obtain the Landau- Ginzburg free energy [lol, the system comes into 
equilibrium on a time scale short compared to the time scales of interest. 

For the remainder of this paper we will take ?7(x, t) and ric{^,t) to be generated by a 
Gaussian distribution with zero mean. That is, <ri{'x,t)> = <i]c{^,t)> = 0, and 

<r/(x, t)ri{x.', t')> = ksTSix - x')S{t - t') (23a) 
<r7e(x, t)r/,(x', t')> = kBTV'^6{^ - x.')6{t - t'). (23b) 



We can use Eqs. (l20l) and (l22l) to determine the time dependence of the decay of fluctu- 
ations in (near-)mean-field systems. The scaling of 0f(x) in Eq. ( ITTj) implies that the cubic 
term in Eqs. fl20l) and fl2T]) can be neglected. A straightforward calculation shows that the 



fluctuations decay exponentially with characteristic times that diverge as e ^ in model A 32 1 
and B?e~^ in model B 32|. For the remainder of this paper we will consider only model A. 



III. PARISI-SOURLAS AND LIFETIME OF FLUCTUATIONS 



The Parisi-Sourlas approach 29|, l30l] begins with the Landau- G inz burg 



equation. Because 



the noise ^^(x, t) in Eq. ( l20l) is Gaussian, the measure of the noise is 29l. l30| 

^^^^ ^ exp[-/3/t/xrftr7^(x,t)] 

/ 5r7exp[— /? / dii.dtrf{y., t)] 

We use Eq. ( 120|) to replace ri{-x,t), let /i = for simplicity, and express Eq. (p4l) as 



(24) 



P{(f)) oc J{(j), rj) exp<—/3 / dx.dt 



^0(x,t) 
dt 



1 2- 



fMA(-i?'VV(x,t)+2e(/.(x,t)+4(/)3(x,t)) , (25) 
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where the Jacobian J(0, 77) of the transformation from rj to (f) is the determinant of the 
operator 6ri{'x,t)/6(f){'x,t). Following Parisi and Sourlas 29|, |30| we introduce the Grassman 
variables ?/'f(x, t) and Tpp^'^.t) which satisfy the algebra 

(26) 

(27) 
(28) 

Because the variables ipp and ipp anticommute, they are referred to as as fermions; 0(x, t) 
is a boson. With this algebra we can evaluate the Jacobian in Eq. (!25i) and write -P(0) as 



'ijjp{x,t) =iljp{-x,t) = jdil)F{y^,t) = jdil)F{y^,t) = 0, 

^ 1pF{y^-it)dll)F{^-it) = 4'Fi'^,t)d'ljjF{^,t) = 1. 



P{(l),1pF,1pl 



exp I - /3 [ Jdxdt (5*5(0, tpF, ^f) + 5*^(0, tpF, ^Pf)] | 



(29) 



where Z is a normalization factor. The quantities Sb and Sf are given by 



SBi(t),'4'F,1pF) 



dyidt 



>90(x,t) 
dt 



+ Ma{- /2'VV(x, t) + 2e0(x, t) + 403(x, t))] ' 



d 



Sf{(P,i1^f,iPf) = I d:s.dt^F{^,t)\^-^ + Ma{ - R^V^ + 2e + 12(j)\x,t))^iPF{^,t). 



(30) 
(31) 



We first consider Sb{4',iPf,'4^f) in Eq- (l30|) . Among the terms found by evaluating the 
square of the term in brackets is the contribution 

a0(x,t) 



C(0) = 2Ma / rfxrft 



dt 



■ [ - i?2vV(x, t) + 2e0(x, t) + 403(x, t 



2Ma / rfx f ^ dt^H{(f){^,t)), 
Jtj 



'dt 



(32a) 
(32b) 



where H is given by Eq. (jll) with h = and (p is replaced by 0(x, t). The integral with 
respect to t gives 

C(0) = 2Ma j d^[H{<p{^,tF)) - H{<f){^,tj))]. (33) 

Parisi and Sourlas assume that tp and tj can be found such that C(0) = and show that 
with this assumption there is a transformation that maps fermions and bosons into each 
other and keeps P{4>) in Eq. (!29l) invariant. They refer to such systems as supersymmetric. 
If a system is in equilibrium, such values of tj and tp can always be found 29|, l30 |. 

The supersymmetric form of P{(f),iljF,'^F) is the proper representation for investigating 
the morphology of the fluctuations in the neighborhood of (near-)mean-field critical points. 
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Near the latter the Hamiltonian in Eq. (jl]) can be assumed to be Gaussian. This assumption 
imphes that the Landau- Ginzburg and Cahn-Hilhard-Cook equations can be hnearized for 
G ^ 1. Because the dependence in the fermionic contribution to the action Sf{4>, V'f, V'-f) 
comes from the nonhnear term in Eq. ( 120|) . hnearization makes the two contributions to the 
action, S'b(0, t/'f, ^f) and 5*^(0, t/'f, V'-f), independent. The integration over the fermionic 
variables '^'^(x, t) and ippi^yt) can be done immediately resulting in the measure 

p^^^ ^ ewj-Pjd^dt i'-^Y + Ml[-i?^VV(x,t) + 60(x,t)]^} ^^^^ 

Zs 



where Zs is the functional integral over (f) of the numerator in Eq. ( ]34l) . 

In equilibrium P(0f) should give the same probability of a fluctuation as PB(0f) in Eq. ([6]). 
To understand the relation between these two probabilities we note that if Pb{4>i) is of order 
e~^, then P{4>{) should also be of order e~^. (This requirement follows from the fact that 
we expect the probability of variations from equilibrium to decay exponentially.) If we take 
(pf in Eq. (IMl) to describe an equilibrium fluctuation, we expect that 



J d:x.dt(^- 



;.f(x,t)x2 



A. (35) 



dt 

where A is a constant independent of 0f. Without loss of generality we can set A = 1. 
Because the spatial extent of 0f scales as the correlation length ^, we have from simple 
scaling arguments that dx in Eq. fl35l) scales as If we use the scaling of 0f(x) in Eq. ([1] 
and the scaling of ^ in Eq. (!9al) . Eq. (l35l) implies 



(36) 



or 

Tf,c ~ e^^- (lifetime of fluctuations) (37) 
Equation fl37p is the well known scaling relation for critical slowing down near mean-field 



critical points for model A (34l |. 



If we require that / di^dt M\ [ - _R^V^0f (x, t) + e0f (x, t)] ^ ~ 1 (see Eq. ([MD), and use the 
scaling relations for ^, (pi, and Tf^c and the same arguments used to obtain Eq. fl36l) . we find 



]^d^2-d/2 

which implies that Ma is a constant of order 1. 



i, (dSj 
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These results are all expected. Note that there is a significant conceptual difference 
between Pb{<P) and P{<f))- The quantity Pb{4') is the fraction of independent members of 
an ensemble in which 0(x) is realized when a measurement is made. Because the system is 
in equilibrium, we can divide a time sequence of measurements into independent segments 
that can be thought of as members of an ensemble. These segments have a duration of 
the order of the decorrelation time (or longer), which is of order Tf ^ near the critical point. 
The quantity P{(f)) gives the probability of a "path" in (p space. The paths of interest here 
are those whose probability is the order of e~^. The path that results in an object with a 
density difference from the background of magnitude 0f = e^^^/G^/^, spatial extent of order 
^ = i?e~^/^, and a lifetime of order is one such path. 

Suppose that in equilibrium there is an object with spatial extent ^ = Re~^^^ but a 
different density. In particular, suppose there is an object of density of 0fc ~ e^/^/G. For 
reasons that we will discuss in Sec. |V]we will call this object a fundamental cluster. Because 
we have assumed equilibrium, we have supersymmetry, and the action is the sum of two 
contributions as in Eq. (1341) . One term has the form as in Eq. (l35l) : 



Ssa=ld.dt('-^)\ (39) 



dt 

We use reasoning similar to that following Eq. fl35p and define the lifetime of the fundamental 
cluster to be given by Sb,i ~ 1- If we substitute 0fc ~ e^/^/G in Eq. (!39l) . we obtain 

or 

Tfc^c ~ 2-d/2 ~ ~q'' (lifetime of fundamental cluster) (41) 

Equation fHTj) gives the lifetime of an object (the fundamental cluster) with density difference 
from the background 0fc ~ e^/^/G near a mean-field critical point. 
The Boltzmann probability Pb{<P) of finding such an object is 

PB(0fc) oc exp^^-pj cix[^[V0fc(x)]2 + e0fc2(x)] }. (42) 

If we use the scaling relations and set /5 = 1 for convenience, we obtain 

PB(0fc) oc e-V^'^^^-'^^l (43) 

We see that Psi't'ic) 7^ P{<Pfc) despite the fact that the system is in equilibrium. If we 
want the probability P{(pic) that there is a path that consists of an object (the fundamental 
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cluster) with density 0fc ~ e^f^ /G, spatial extent ^ ~ -Re and lifetime e ^ rather than 
e'^/G, then using Eq. (l34l) we have 

P(0fc) ~ e-^/^^^^-^^^ (44) 

and PB{<Pic) = P{4'ic)- The implication of these results is that the Boltzmann probability 
Pb{'Pic) requires a given time, that is, the decorrelation time, which near a mean- field critical 
point scales as e^^. The probability P{4>ic) = Pb{'Pic) only if t is chosen to be the decorre- 
lation time Tf c- In general, the decorrelation time is not equal to the lifetime of the object 
of interest. Note that the same arguments apply to the normalization factors Z in Eq. ( fT2|l 
and Zs in Eq. In particular, Z = Zs only if the time scale is chosen to be r ~ e^^. 

The mobility Ma need not be a constant independent of e [3]. For an object with density 
(/)fc ~ e^/^/G, the second term in the action in Eq. fl25l) has the form 

^2 = -j dxdt Ma^ [i?2[V0fc(x, t)]2 + 2e0fc(x, t)] \ (45) 



If we use the scaling relations and the lifetime given by Eq. fl4T|) . we obtain 

MA^^R'^e-^/h-^ 



1, (46) 



or 

Ma ~ i?^e2-^/2 > 1. (47) 

Ma in Eq. (H7I) depends on e in contrast to the mobility in Eq. (1381) . If we had considered 
a lifetime of rather than e^^/G in Eq. ( |4T|) . would be order unity. 

In summary, the probabilities Pb{4>) and P{(f)) are equal if the lifetime of an object is 
the order of e^^ near the mean-field critical point. The lifetime of an object is obtained by 
requiring that P(0) ~ e~^. For objects with density e^^^/G and a lifetime of e~^/G, the 
usual Boltzmann factor will not give the probability of observing such an object. 



IV. SPINODALS AND PSEUDOSPINODALS 



In this section we discuss the meaning of spinodals and pseudospinodals. We begin 
with the Hamiltonian in Eq. (jl]) and the partition function in Eq. ([5]). For G ^ 1 the 
partition function can be evaluated using saddle point techniques, and the free energy has 
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the Landau- Ginzburg form 28| . 

F=Jd^ J:|![V0(x)]2 + e02(x) + 0'*(x) - /i0(x)] . (48) 
We set the gradient term equal to zero to obtain the free energy density 

For e > there is only one real extremum of the free energy. For e < there are three real 
extrema, one maximum and two minima. For h = there are two states or values of with 
the same free energy. As \h\ is increased, one of the minima becomes higher than the other. 
The higher minimum corresponds to the metastable state. Increasing \h\ further eventually 
results in the disappearance of the metastable minimum. This value of \h\ is referred to as 
the spinodal field hg. It is easy to see from Eq. fj49|) that a.t h = hg, f has an inflection point 
at = (ps- If we set (for e < 0) df /d(j) = d'^f/dcj? = 0, we find 

0, = (|e|/6)i/2 and = 4|e|=^/V(3v^). (50) 

We define the new variable Ah, 

Ah = hg- h, (51) 

and the new field, ipi'^) = — 4>s + cl, and write the mean-field free energy as 

F=J dx[^[V^(x)]2 + Ah'/^X^^\^) - A2^=^(x) + A3^^(x)]. (52) 

The parameter a is chosen so that the term linear in '0(x) does not appear in F. The 
coefficients Aj are functions of e and independent of Ah. (More precisely, the Aj are a 
function of Ah, but as Ah 0, the Aj approach constants.) The free energy in Eq. (!52ll is 
constructed so that the spinodal is at -0 = and Ah = 0. 

As for the critical point we assume that the fluctuations associated with the spinodal 
can be described by a Gaussian-Landau-Ginzburg- Wilson Hamiltonian with the partition 
function 

Z=j6^exp [- [3 J d^^[V^{^)f + Ah^/^Xi^^i^)]. (53) 
If we follow the same argument that we used at the critical point in Sec. [Tll we obtain 

C ~ RAh-^'^ (54a) 
X ~ Ah^^'^. (54b) 
The fluctuations of the order parameter density scale as 
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where the Ginzburg parameter near the spinodal is 

Gs = R'^h"^-'''\ (55) 

The system is mean-field when Gs oo and near-mean-field for ^ 1. A spinodal is 
only a true critical point in the mean-field limit . The distinction between a spinodal 

{Gs oo) and a pseudospinodal {Gg ^ 1) will be clear from the context. 

As discussed in Sec. [Tll the scaling of the order parameter density and the density of the 
order parameter fiuctuations is not the same in (near-)mean-field systems. For example, 
near the critical point for e < we have from Eq. fH9|) 

-2|e|0 + 403 = 0. (56) 

Equation (l56l) gives the order parameter density at the minima for e < 0. As |e| — > we 
have 

0~|e|i/l (57) 
Near the spinodal we use Eq. (132]) and assume that ipi^x) is independent of x. We have 

f{,p) = Ah'/'Xiij' - A2^' + Aa^^ (58) 

As Ah —>■ 0, scales as 

-0 ~ Ah^^"^ (order parameter density near the spinodal), (59) 

where we have dropped the ip'^ term in Eq. (1581) . Equation (IS^ gives the scaling of the order 
parameter density. As near the critical point, the scaling of the order parameter density and 
the density of the fiuctuations is not the same in mean-field systems near the spinodal. 

In the above discussion we kept the temperature fixed and approached the spinodal 
by varying the magnetic field h. Alternatively, we can keep the magnetic field fixed and 
approach the spinodal by varying the temperature. To obtain the critical exponents in the 
temperature variable we return to Eq. ( 149|) and write e as 



T-T^ T-Ts Ts-Tc 

^ ^ = ^ + ' =e + As, (60) 
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where is the spinodal temperature for a fixed field h = hg. We write 



= -26 -2\As\ + 12(j)i + 24(f)s^ = 0, (61) 

's 

Because and (ps are on the spinodal curve, we have 

-2|A,| + 1202 = 0, (62) 

-2|£| + 240,V^ = 0. (63) 

From Eqs. (15^ and we have 

~ £, (64) 

which implies from Eq. fl54cp that as the spinodal is approached, the density of the order 
parameter fiuctuations scales as 

M^^O- [^,g3-d/2]i/2 - (65a) 
Similarly, the correlation length scales as 

e ~ Re-^^\ (65b) 

and the suceptibility diverges as 

X ~ (65c) 

Equations and fl65|) give the critical exponents near the spinodal in terms of Ah and e. 

We first discuss the application of the Parisi-Sourlas method near the spinodal. If we 
construct a Landau- Ginzburg equation from the free energy in Eq. (!52l) . we have 



-Ma,s [ - i?'VV(x, t) + 2A/i^/2^i^(x, t) - 3A2^'(x, t) + 4A3V^^(x, t)] + r/(x, t). 



g^(x,t) 

dt 

(66) 

Because the noise 77(x, t) is Gaussian, we obtain an expression for the probability of a path 
?/'(x, t) of the form 

exp [- Pjd^dt {'-^Y + Ml{ - R^V^jji^, t) + 2A,A/^V2^(x, t) }^] 
^sp(^/') = — (67) 

where we have used the linear form of the Landau- Ginzburg equation and have assumed that 
the system is in metastable equilibrium, which implies supersymmetry. Arguments similar 
to those used at the critical point show that the relaxation or decorrelation time scales as 

r ~ A/i^^/2. (decorrelation time near the spinodal) (68) 
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r.y,. 



FIG. 1: The inverse susceptibility as a function of the magnetic field h (from Ref. lZ|). Note that as 
the number of neighbors q is increased, the inverse susceptibility more closely follows a power law. 
The inset shows the behavior of (ksTx)^^ closer to the pseudospinodal. 



All considerations of the difference between equilibrium Boltzmann probabilities and prob- 
abilities of paths are the same near spinodals and critical points. 

We now consider the nature of the pseudospinodal. As mentioned, for ^ 1 but 
finite, there is no spinodal. However, the system behaves as if it existed if Gs is sufficiently 
large. In Fig. [1] we plot the inverse of the isothermal susceptibility x found by a Monte 
Carlo simulation for a d = 3 Ising model as a function of the applied magnetic field h for 
different values of i? [t]. The temperature is taken to be 4Tc/9, where Tc is the critical point 
temperature. The solid line is the mean-field prediction for q — > oo, where q is the number 
of spins that interact with a given spin 3^. Data was taken only if the metastable state 
lived longer than 10^ Monte Carlo time steps per spin. For nearest-neighbor interactions 
{q = 6) the data stops at /i ~ 0.5 far from the spinodal value of the field hg = 1.43. As 
R and hence q is increased, the data approaches the mean-field result and the spinodal can 
be more closely approached. This result indicates that the larger the value of R, the more 
system behaves Uke there is an underlying spinodal. 

Another way to understand the nature of the pseudospinodal is to look at the behavior 
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of the zeros of the partition function as a function of R. The zeros of the partition function 

corresponding to the spinodal he in the four- dimensional complex magnetic field-temperature 
P 

space for finite R [9]. As R is increased, the zeros move toward the real {h,T) plane similar 
to the behavior of the zeros of the partition function for Ising models in finite systems as 
the system size increases jsol, [s^]. The idea is that the pseudospinodal appears to be a 
critical point if h is not too close to hg. What is meant by too close can be estimated by the 
magnitude of the Ginzburg parameter Gg in Eq. fl55|) . The value of Ah where the spinodal 
concept fails can be made smaller by increasing R. Hence, the theoretical arguments we 
made about the properties of fluctuations near the spinodal can be tested in systems where 
the interaction range R is large, even though there is no true spinodal in nature. However, 
such statements have to be modified in systems with a phase transition that involve spatial 
symmetry breaking such as the liquid-solid transition (see Sec. IVIip . 



V. FLUCTUATION STRUCTURE 

In this section we will use scaling arguments and the cluster mapping discussed in Ap- 
pendix |A] to determine the structure of the fluctuations in (near-) mean-field systems. 

Ising critical points and the spinodal in (near-)mean-field systems have been mapped onto 



percolation transitions 25|, 138|. For simplicity, we will assume that the interaction between 
spins in is a constant up to a distance R and is zero for distances greater than R. To map 
the critical point onto a percolation transition we toss bonds randomly between pairs of 
parallel spins that are separated by a distance less than or equal to R with a probability 
Pb = 1 — e~^^'^ where J is the usual Ising coupling constant. This mapping guarantees that 
the percolation transition occurs at the critical point. For the spinodal the bond probability 
is = 1 — e~^'^^^~P\ where p is the density of spins in the stable state direction; that is, 
if the metastable state is in the up direction, then p is the density of spins in the down 
direction. The size of the clusters is determined by the number of spins in a cluster. The 
details of this mapping are given in Appendix |X1 

Equilibrium critical phenomena fluctuations in the order parameter are defined as devia- 
tions from its mean value with a linear dimension of the correlation length and a free energy 
cost of order 1. These properties were used in Sec. II to calculate the density of fluctuations 
near a critical point. We will see that the same properties apply near a spinodal. 
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The geometrical quantity that is isomorphic to the free energy is —ksT times the mean 
number of clusters (see Appendix |A]) . The free energy density near the critical point scales 



as e^, because the specific heat exponent a = for the mean-field critical point 32|, l33 |. 
Hence the free energy F{e) in a correlation length volume scales as 

F(e) ~ e^e'^ = R''^-'''\ (69) 

and the mean number of clusters in a correlation length volume scales as 

rlfc,c ~ R^e^-'^'^ = G. (70) 

How are the clusters related to the fluctuations? In non-mean-field systems such as the 
Ising model with i? ~ 1 and d < 4, the mean number of clusters in a correlation length 
volume near the critical point scales as e^-a^d ^ ]^ ^^le isomorphism between the Ising 
model and percolation implies that the the pair distribution function p^"^^ is the same as the 

(2) 

pair connectedness function pc , which is the probability that two spins a distance x apart 
belong to the same cluster. For x ^ ^, pc^^ is roughly a constant and equal to ^^j, where 
0ci is the density of spins in the cluster. This density must be equal to the density of a 
fluctuation which scales as (see Eq. fl9bl) ). Hence, 0ci scales as in a non-mean-field 



system, and the clusters are a statistical realization of the fluctuations |38l. |39|. 

For G ^ 1 the number of clusters in a correlation length volume near the critical point 
scales as G (see Eq. fITOl) ). which is much larger than unity. Are the clusters a statistical 
realization of the fluctuations? To understand that the answer is no, we note from Eq. (fT7|) 
that the density of a critical phenomena fluctuation scales as 0f(x) ~ e^/^/(i?'^e^~°'/^)^/^. If 
a single cluster were to correspond to a fluctuation, then the density of spins in a correlation 
length volume would scale as 

G0,(x) = R'e'-'^'j^j^^, = [R^e-'^r^e^'^- (71) 

Because the mean-field limit corresponds to letting R oo before e — > 5|], Eq. fl7T|) 
implies that the spin density is infinite in the mean-field limit, which is impossible. Hence, 
the density of the clusters must be much smaller than the density of the fluctuations, and 
a fluctuation does not correspond to a single cluster as it does in short-range systems. We 
conclude that the clusters in (near-)mean-field systems play a different role, and we will refer 
to the clusters in these systems as fundamental dusters. 
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To understand the relation between the fundamental clusters and the fluctuations in a 
(near-)mean-field system, we again use the the fact that the pair distribution function p^"^^ 
is isomorphic to the pair connectedness function pc . Because the latter is the probability 
that two sites a distance x apart belong to the same cluster, we have 

pf(x<0-Pfc.c^||^^, (72) 

where pic,c is the probability that the first site belongs to any one of the i?'^e2-d/2 d^gters, 
and pfc,c/-R'^e^~'^^^ is equal to the cluster density, which is the probability that another site 
belongs to the same cluster as the first. 
Similarly, we have that 

.1/2 

(^dg2~2)l/2 

We have pf\^ ^ = P'c\^ ~ 0; ^.nd hence We have 



,2 



Pic 



c 



.1/2 



1 2 



^Rd^2-d/2^ 



1/2 



(74) 



and pfcc = e^^^- Hence, the density of spins in a fundamental cluster is 



Pic c e^/^ gi/2 

0fc,c(x < ~ j^dJ-d/2 = Rd^2-d/2 = (75) 

Because the density of the fundamental clusters is much smaller than the density of the 
fluctuations for G ^ 1, a fluctuation must correspond to many fundamental clusters. 

We can test the prediction for 0fc,c(x) in Eq. fl75|) by determining the dependence of mfc,c, 
the mean number of spins in a fundamental cluster, on e. This dependence is given by 



In Fig. [2] we plot nid a function of e for fixed R. The slope of the log-log plot is consistent 
with the theoretical prediction in Eq. fl76l) . 

We now discuss the relation between the fundamental clusters and the fluctuations in 
more detail. As discussed in Appendix |A] the clusters are constructed to be independent. 
Therefore a given cluster can "flip" independently of the other clusters. There are i^'^e2-d/2 
clusters near the mean-field critical point, half up and half down by symmetry. Because 
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FIG. 2: The mean number of spins in a fundamental cluster with spatial extent ^ found in a 
Monte Carlo simulation of the d = 2 Ising model with R = 20 near Tc as a function of e. The linear 
dimension of the system is L = 240. The slope is ~ —1.5 if the data is fitted in the range [0.03, 0.2]. 
Note the deviation of rrid from the mean-field prediction near the mean-field critical point where 
G is too small to apply mean- field arguments and where finite size effects become important. 

the clusters are independent, the mean number of excess fundamental clusters in a given 
direction is determined by a random walk: 



^fc,c 



(77) 



More precisely, the distribution of the number of fundamental clusters in a fluctuation is a 
Gaussian peaked about Afc,c- 

From this analysis we see that the density of a fluctuation is the product of the density 
of fundamental clusters, Eq. flTSl) . and the mean number of excess clusters, Eq. ( 1771) : 



.1/2 



<0 ~</>fc,cAfc,c ~ j^^^2-d/2 



.d 2-d/2\l/2 



.1/2 



a/2 



(78) 



' (i?dg2-d/2)l/2 ^1/2' 

in agreement with Eq. 0171) . 

A similar analysis can be done near the spinodal, and we will only summarize the results 
here. In this case there is an infinite cluster that is a statistical realization of the metastable 



state magnetization 23|, |2J]. If we subtract this cluster 40|, the results are similar to those 
near the mean-field critical point. The mean number of fundamental clusters in a correlation 
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length volume, half up and half down, is given by (compare to Eq. (1701) ) 
and their density is (see Eq. fl75l) ) 



(79) 



(80) 



RdAh^/2-d/4 ■ 

We see that a system is near-mean-field when the number of fundamental clusters in a 
correlation length volume is large (see Eq. ( !79l) ). As e{Ah) is decreased for d < 4 (critical 
point) or d < 6 (spinodal) for fixed R, G {Gs) decreases and the system becomes less mean- 
field. 

If we use the same random walk argument as for Eq. fl77|) , we find that the density of the 
order parameter fluctuations near the spinodal scales as 



< 







(81) 



(i?^A/l3/2-d/4)l/2 ^1/2 • 

Hyperscaling (two exponent scaling) is satisfied for systems that are non-mean-field if e 



(A/i) is sufficiently small 32, ISSj. This connection suggests that the magnitude of G{Gs) 
determines the existence of hyperscaling. If G {Gs) decreases, the clusters must coalesce 



and additional length scales are introduced 



4l|. We now show that the existence of one 



relevant or divergent length scale is insufficient for the existence of hyperscaling if G ^ 1 
because the number of fundamental clusters changes as the critical point is approached. The 
assumption of one divergent length scale leads to the following form for the singular part of 
the free energy density near the critical point 32, 



If we differentiate / twice with respect to h and set = 0, we find 

a2/(e, h) 



h=0 



(82) 



(83) 



The left-hand side of Eq. ( l83l) is the isothermal susceptibility x- We now fix ^^^e to be 
equal to one. Because /(I, /i = 0) is not singular 32l, l33| and ^ ~ g-i/s/T^ have ^'^y^-'^ = 



-(2y,,--d)/yj 



. Hence the exponent 7 that characterizes the divergence of x near the critical 



point is given by 



7 



2yfe - d 



(84) 



By using a similar argument, we obtain 
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/3 = (85) 
Vt 



For mean- field Ising models or simple fluids (3 = 1/2 and 7 



32 



. Hence i/h = 3d /A and 



33I ] and not d/2. Hence, 



Vt = d/2. For fixed R, ^ = Re'^'^, so that yr = 2 for all d 
for fixed -R two exponent scaling does not hold in the neighborhood of a mean-field critical 
point. The same argument holds near the spinodal. 

In contrast, consider what happens for fixed G = R'^e^^'^^^. From Eqs. (|T7j) and ( JTSll we 
have 

r ,^V2 -1 2 

If we keep G = /^"'g2-a!/2 constant, then R^e"'^^'^ oc e~^. Therefore the susceptibility x 
so 7 = 1. Likewise from Eq. (1751) the density scales as e^/^ so P = 1/2. For G fixed 
i?e~^/^ oc e"^/'^ so that v = 2/d and ?/t = (i/2. Hyperscaling now holds and 7 + 2/5 = (ii/. 
Similar arguments hold near the spinodal if the infinite cluster is removed. 

In summary, we have shown that the relation between the clusters and the fluctuations 
in (near-)mean-field systems is more complex than in non-mean-field systems. In particular, 
the individual clusters are not realizations of the fluctuations, which instead are related to 
fluctuations of the number of clusters. For this reason we refer to the clusters in (near- 
)mean-field systems as fundamental clusters. The mean number of fundamental clusters in 
a correlation length volume is proportional to the Ginzburg parameter G. The dependence 
of G on e {Gg on A/i) causes the breakdown of hyperscaling. 



VI. LIFETIME OF THE FUNDAMENTAL CLUSTERS 

In Sec. |V] we assumed a scaling form for the free energy. Here we do a more detailed 
calculation using the Landau-Ginzburg- Wilson Hamiltonian in the mean-field limit. Our 
main result is the scaling dependence of the lifetime of the fundamental clusters. We also 
recover the same scaling results for the free energy. 

If we set h = 0, scale all lengths with the correlation length, and assume that 0(x) scales 
as e^/^ near the critical point, we obtain the Hamiltonian in Eq. ([8]). If we assume that 
0(y) is independent of y and restrict the integral to a region the size of a correlation length 
volume, we have 

= i?V-'^/2[±02^04]^ (87) 
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and the partition function becomes 

Z{e)=l 4e-^«'^'-'''(±^'+^'). (88) 



For G = R^e^-'^'^ > 1 we can do the integral in Eq. (IHHl) using saddle point techniques. For 
e < the saddle points are at = ±l/-\/2, and we obtain 

Z(6) oc e^«'^"''^ (89) 

and hence the free energy is 

- kBT\nZ{e) = -R^t^-^'^. {G > 1) (90) 

We have neglected the logarithmic corrections generated by the steepest descent integral. 
Note the minus sign on the right-hand side of Eq. ( |90l) . For e > the saddle points are at 
±.i/\/2 and the free energy is also proportional to — i?'^e^~'^/^. Hence as argued in Sec. IVland 
Appendix |X] we can use the percolation mapping to show that the number of fundamental 
clusters scales as i?'^e2-d/2 ^^^^ ^^iq critical point. As usual, similar arguments can be used 
near the spinodal for — > oo and near the pseudospinodal for ^ 1. 

To determine the lifetime of the fundamental clusters, we return to the Parisi-Sourlas 
method. At the spinodal we have from Eq. ( 1671) 

2 



1. (91) 
As in Sec. IIIII the lifetime of a fundamental cluster is found by requiring that 

(92) 



or 

Tfc s ~ T^TT-TiTTj— 777- (fundamental cluster lifetime near spinodal) (93) 

We see that near the critical point (G ^ 1) and the pseudospinodal {Gs ^ 1), the lifetime 
of the fundamental clusters is considerably shorter than the lifetime (decorrelation time) of 
a fluctuation near the critical point (e~^) and the spinodal (A/i~^/^). 

To understand the relation between the lifetime of the fundamental clusters and the 
lifetime of a critical phenomena fluctuation recall that the clusters are independent. We 
consider the fluctuations to be formed from the "vacuum" (zero magnetization near the 
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critical point and zero net magnetization after the infinite cluster is subtracted near the 
spinodal) by a random walk in the number of fundamental clusters. At the critical point the 
density of the critical phenomena fluctuations is given by Eq. ( ITTl) . Because the fluctuations 
arise from a random walk in the number of fundamental clusters, there must be a "walk" of 
(i^'^e^"'^/^)^/^ cluster flips (steps) in the direction of the fluctuation to obtain a density of 

pl/2 

The time needed for ]^d^2-d/2 attempted cluster flips is 

in agreement with Eq. fl57|) . The same considerations near the spinodal yield Tf s ~ A/i^^/^ 
for the lifetime of a fluctuation, in agreement with our earlier result for Tf.s in Eq. fl68l) . 

In summary, we have argued that the fluctuations near a mean-field critical point and a 
spinodal are not represented by a single fundamental cluster. The relation between these 
clusters and the fluctuations in (near-) mean-field systems is qualitatively different than in 
systems that obey hyperscaling. In the former the fluctuations are formed by a random 
walk in the number of fundamental clusters that "flip" on a time scale much shorter than 
the scale set by critical slowing down. A summary of our notation and our main results so 
far is given in Table [H 

VII. CLUSTER STRUCTURE AND INSTABILITIES IN SUPERCOOLED 

LIQUIDS 

We begin our discussion of the consequences of the fluctuation structure in (near-)mean- 
field systems by considering the liquid-solid spinodal in supercooled fluids. To explain the 



role o 
wood 



the structure of the fluctuations we first provide some background. In 1951 Kirk- 
noted that approximate equations for the distribution functions in the liquid state 



appeared to show an instability as the supercooled liquid is 



quenched deeper. Kirkwood 



began with the first equation of the static BBGKY hierarchy |43l] 

- kBTVip'-'^^i) = j rfx2 Viy(xi2)p(')(xi,X2), (96) 

where ^'•^-'(xi) and p*^^''(xi, X2) are the one and two particle distribution functions respec- 
tively, Vi denotes differentiation with respect to the position of particle 1, and the interaction 
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quantity 


mean-field critical point 


spinodal 


hyperscaling 


Ginzburg parameter 


G = Rde^-dn 


Gs = i?"'A/i3/2-'^/4 


fixed 


order parameter 








fluctuations in order pa- 
rameter density 


<Af (x < ~ eVVGi/2 


V'/(x<e)~ei/VGy' 




lifetime of fluctuation 








lifetime of fundamental 
cluster 




''"fc.s ^ /Gs 




mean number of clusters in 
correlation length volume 


nfc,c ~ G 


^fc,s ~ Gs 


~ 1 


density of fundamental 
cluster 


0fc,c(x < e) ~ e'/'/G 


V'fc,s(x < e) ~ e^^VG, 




mean number of excess 
fundamental clusters 




Afc,s ~ G^^ 





TABLE I: Summary of our notation and some of the important scaling relations derived in the text. 
The spinodal can be approached by reducing the magnetic field difference A/i for fixed temperature 
or by decreasing the temperature difference e for fixed magnetic field. The two approaches are 
related by A/i^/^ ~ e. The exponent z characterizes critical slowing down and is of order 2 for 



systems described by model A 



34|. 



potential V{xi2) is assumed to be pairwise additive and spherically symmetric. Suppose that 
the system is in the Uquid phase where p'^^^(xi) is a constant equal to p, and p(^''(xi,X2) is 
a function of Xi2 = |xi — X2_| and is equal to p^h^xu) = p^(l + g{xi2)), where g{xi2) is the 



pair correlation function [43|. We substitute 

p(^)(xi) =p + c<;(xi) 

into Eq. ( l96|l . treat a;(xi) as a small perturbation, and linearize Eq. ( l96l) to find 

,Viu;(xi) 



-ksT 



P 



(iX2 VlV{Xi2)h{xi2)uj{^2), 



where the spherical symmetry of V{xi2) and h{xi2) results in 

dx2ViV{xi2)h{xi2) = 0. 



(97) 



(98) 



(99) 
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We have also assumed that a possible instability in h{xi2) is higher order 'm. uo. If we define 

g(xi - X2) = ViV{xi2)h{xi2). (100) 



we see that there is an instability 42| if there is a nonzero solution to 



fcBTViCo'(Xi) + Pp rfX2 g(Xi - X2)co'(x2) = 0. 



(101) 



Kirkwood analyzed Eq. (110 II) for the hard sphere fiuid and found that there was an instability 



in (i = 3 42|. However, he ignored a possible instability in h{xi2) which is related by the 
BBGKY 43| hierarchy to possible instabilities in all of the distribution functions. Hence, it 
is not clear that the instability Kirkwood found is real. A more careful analysis [4^ suggests 
that the instability vanishes when higher order terms are considered. 

To investigate the existence of an instability and its relation to a possible spinodal, Grewe 



and Klein 



45 



46l ] investigated the properties of a simple fluid for which the interaction 



potential has the Kac form ^] given in Eq. ([T]) with Vr = and 



$(7|x|) 





if 


7 X 


< 1 




if 


7 X 


> 1 



(102) 



45 



46| 



In the mean-field limit 7 0, it was shown that all distribution functions of order higher 
than two are completely specified by only the single particle a nd p air distribution functions 
and that p'^^)(xi) in the limit 7 — > satisfies the equation 

p«(xi) = zexp[-P j rfx2$(|xi2|)p«(x2)], 

where z = e~^^ and p is the chemical potential. Similarly, (7(|xi2|) satisfies 

^(|X12|) = /3p$(|xi2|) -Ppj rfX3^(|xi - X3|)<I>(|X2 - X3I), 



(103) 



(104) 



where all length scales are in units of 7"^ = R. Note that o(|xi2|) is of order 7"^. The 
derivations of Eqs. (11031) and (11041) are given in Refs. [45! and |46. 

From Eq. (I104p the structure function S{k), which is obtained by taking the Fourier 
transform of (7(1x12!), is proportional to 

1 



S{k) oc 



(105) 
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where $(fc) is the Fourier transform of $(7|x|) and k = |k|. Note that the structure function 
is order one in the 7 — > hmit. 

We can perform a stabihty analysis on Eq. fllOSp similar to that done by Kirkwood on 
Eq. fl96|) . We substitute Eq. fl97|) into Eq. fllOSp and linearize in ci;(xi) to obtain 



UJ 



xi) = -Ppj $(|xi - X2|)u;(x2), (106) 



where p is the solution of p = 2;exp[— /3<l>(0)p], and $(0) = / (ix$(|x|) > 0. 

There is an instability only if there is a nonzero solution to Eq. (11061) . If we take the 
Fourier transform of Eq. (11061) . we can express the instability condition as 

l + /?pl'(|k|)<0, (107) 

or $(|k|) < for some value of |k|. This condition is satisfied for the potential in Eq. (11021) . 
Because $(7|x|) in Eq. (I102p has a Fourier transform that is bounded from below, there is 
a value of (5p below which there is no instability. If = |k| is the location of the global 
minimum of $(|k|) < 0, then the system has no instability for jSp < — l/<l'(A;o). We see from 
Eq. (llOSp that the structure function S^ko) first diverges for fixed p (as T is decreased) at 
the same value of the temperature T at which an instability first appears. The divergence 
of S{ko) implies that the instability in the mean-field system is a spinodal and is analogous 
to the divergence of the susceptibility at the Ising spinodal. 

In the mean-field limit no higher order distribution functions need to be considered and 



the results of Grewe and Klein 



45 



46| are rigorous. The structure function S{kQ) diverges as 
(T — Ts)~^ so that the critical exponent is the same as the Ising spinodal if the temperature 
rather than the magnetic field is used to approach the spinodal in the Ising model. The only 
difference is that S{k) diverges at = in Ising models and gases rather than at ko 7^ 0. 



45 



46| 



The other critical exponents are also the same as for the Ising spinodal 

We next discuss an important difference between measurements of the spinodal exponents 
in Ising models and in supercooled fluids. As we have discussed (see Fig. [1]), there is no 
spinodal in an Ising model for R finite, but we see spinodal-like behavior for i? ^ 1 if the 
system is not quenched too deeply into the metastable state The larger R, the 

more the pseudospinodal behaves like a true spinodal. 

In the supercooled liquid there is no direct evidence of a spinodal or pseudospinodal from 
either experiments or simulations. We will see that this lack of direct evidence is due to 
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the structure of the fluctuations in (near-) mean- field systems and the crucial role of the 
fundamental clusters. 

For the potential in Eq. fll02p the system is a fluid for high temperatures and/or low 



densities 



45 



46 



471]. If the temperature T is lowered at a fixed density, the liquid-solid 



instability is encountered at the spinodal temperature Tg. If T is lowered below Tg, the 
uniform density fluid phase becomes unstable and a "clump" phase is formed 47|. Because 
we are interested in the nature of the spinodal, the behavior of the system for T < is not 
of interest here. 

Unlike Ising/Potts models, there is no precise definition of a cluster in a continuum 
system of particle. However, it is reasonable to assume that the scaling behavior of the 
fluctuations and fundamental clusters near the liquid-solid spinodal is the same as near the 
Ising spinodal. This assumption is consistent with the fact that the Ising and spinodal 
exponents for the system defined by Eq. fll02p are the same. For convenience, we will use 
temperature scaling near the liquid-solid transition rather than the analog of magnetic field 
scaling. From Eqs. (|93l) and (1641) the lifetime of the fundamental clusters as £ — > scales as 



-1 



^f^''^ ~ j^d 3-d/2 • (lifetime of fundamental cluster near spinodal) (108) 



Grewe and Klein 



45 



46l | used the mean- field formalism developed by Kac et al. In 



particular, the interaction range R = is taken to infinity before the spinodal is ap- 
proached, that is, before the limit e oc (T — Tg) goes to zero. Hence, in the mean-field limit 
there exists fundamental clusters with probability of order one but zero lifetime (see the 
discussion in Sec. IIII|) . In experiments and simulations, the converse is true. For exam- 
ple, consider a measurement of the structure function S{k) in a simulation. The structure 
function is obtained by computing 

j 

where Xj is the instantaneous position of particle j and (■ ■ ■ ) denotes an ensemble average. 
Because a simulation can be performed only on systems with finite R, the time of the 
measurement, which is instantaneous, is much less than the nonzero fundamental cluster 
lifetime in Eq. (11081) . Therefore the result of the simulations need not be consistent with 



the mean- field predictions 45|, |46|. (The time scale of the measurement does not refer to 



28 



the time over which data is taken, but to the time during which the probe (say a neutron) 
is in contact with a fluctuation.) 

Before we consider the behavior of S{ko) in supercooled liquids, we discuss the measure- 
ment of S{k) near the Ising spinodaL We will find that the measured behavior of S{k = 0) 
agrees with the mean-field predictions [ssl. The application of the scaling argument in Sec.lVl 
to the susceptibility near the Ising spinodal gives (see Eqs. (ITSl) and (IHTl) ) 



L(i?^A/l3/2-d/4)l/2j 

Alternatively, we can calculate x directly from the fundamental clusters. The density of 
a fundamental cluster is given in Eq. flHOj) . Hence, the isothermal susceptibility associated 
with one fundamental cluster is (see Eq. (IHTIl ) 



R Ah' ' ~ fill) 



Because the clusters are independent, the isothermal susceptibility of the system is the sum 
of the individual cluster isothermal susceptibilities. Because there are 

clusters 

in a correlation length volume, x = A/i^^/^, consistent with Eq. (11 101) . Although there 
are clusters of both up and down spins, each cluster has the same isothermal susceptibility 
because the susceptibility is proportional to the square of the spin density 
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Another way of obtaining the same result as in Eq. f lllOp is to calculate the structure 
function S'ifc(k) for one fundamental cluster from the Fourier transform of the connectedness 
function of the cluster. We can write (see Eq. (IHOj) ) 



A/ii/2 



2 



^"c,s(k) - [ ^,^^3/2-./4 j (112) 

The delta function comes from integrating over the infinite size of the spatially uniform 
cluster. For a cluster with the spatial extent of the correlation length S,-, 5(k) in Eq. flll2p 
would be replaced by a function whose height is ^'^ and whose width is proportional to 

The difference between the fundamental cluster lifetime rfc,s in Eq- fl93l) and the lifetime 
of a fluctuation ~ A/i~^/^ suggest the following picture. A fluctuation is a collection of 
fundamental clusters that appear and disappear on the time scale rfc,s- The lifetime of a 
fluctuation is a much longer than a fundamental cluster for ^ 1, which implies that the 
fundamental clusters come and go with different angular orientations but with their centers 
in roughly the same place for a time of order Ah'^^^, the lifetime of a fluctuation. If we 
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make a measurement on a time scale ^ A/i~^/^, the number of up and down fundamental 
clusters would be the same on the average, and the measured structure function would show 
no correlations between the spins. 

Consider a measurement on a time scale tmeas such that tmcas ^ Tfc,s and tmcas ~ A/i^^/^. 
Because the measurement time is comparable to the lifetime of the fluctuation, the mea- 
surement will see a spin density equal to the fluctuation density, which is generated by 
fluctuations in the number of fundamental clusters. That is, the individual fundamental 
clusters cannot be distinguished, and 

r Ah^/^ 12 
^(k) - [ ^R.^hS/2-,/.y/2 \ ^(k)- (113) 

If we replace 6(h) by we obtain the same result as in Eq. (11 101) . obtained by considering 
the fluctuations directly. 

Now suppose that a measurement is made such that tmcas ^ Tfc,s- In this case an external 
probe or a simulation would see a set of order R'^Ah^^'^~'^^^ frozen fundamental clusters. 
To determine the structure function that would be measured, we need to add the cluster 
structure function in Eq. (11121) for the R'^Ah^^'^~'^^^ frozen clusters. To do so we convert 
the sum to an integral by using one of the factors of 1/ R^Ah^'^-'^/^ in Eq. (fTT2|l to create 
an infinitesimal element of solid angle dVL. (We ignore numerical factors because we are 
interested only in the scaling properties.) In so doing we are assuming that the R^Ah^/'^-'^'^ 
fundamental clusters overlap each other with random orientations. Hence the sum over 
clusters becomes an integral over solid angle: 

(114) 

For = 3 we have in spherical polar coordinates 

(115) 



Hence 



Sik) - ^3A/,3/2-3/4 ^- (116) 



If we now replace 5{k) by ^ with A; ~ ^ ^, we obtain 
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in agreement with Eq. flllOp . The generahzation to arbitrary dimensions is straightforward. 
We conclude that we obtain the same scahng dependence for the susceptibihty (equal to 
S{k = 0)), independent of the relative order of magnitude of the measurement time and 
consistent with simulations of the Ising model 

We now discuss the measurement of S{k) near the liquid-solid spinodal. As stated, we 
assume that the fundamental clusters in the supercooled liquid scale the same way as they do 



in Ising models. Because the clusters are independent [38|], the structure function of a single 
fundamental cluster has to contain information about the symmetry of the instability. That 
is, a collection of independent clusters cannot generate a symmetry that does not already 
exist in each cluster. Because we are interested in the limit of stability of the supercooled 
liquid and know that the instability occurs at ko ^ 0, the clusters must reflect this symmetry. 
Hence, we will assume that the clusters have a symmetry reflected by the wave vector kg 
with arbitrary orientation. In analogy with the Ising spinodal, we expect that the structure 
function S'fc,s(k) of the fundamental clusters can be approximated near the spinodal by 

^^c.(k)~[^5ifc^]V-ko). (118) 

We have used temperature variables rather than the chemical potential, the analog of the 
magnetic field. There are other peaks in 5'fc,s(k) at |k| 7^ ko, but we will focus on 5'fc,s(ko), 
the peak associated with the divergence as the liquid-solid spinodal is approached. 

We first consider a measurement on a time scale such that tmeas ^ Tfc.s, where rfc^s is 
give by Eq. fllOSp . As before we need to sum over all orientations of the frozen clusters 
whose centers are fixed. We convert this sum to an integral by absorbing one of the factors 
Eq. (IllSp to form an infinitesimal. In this case we need to integrate over 
orientations of the vector ko keeping the magnitude ko constant. For c? = 3 we have 



-2 



1 



^(^) ~ ^3^3-3/2 J ded<p6{k - ko)6{e)6{<P). (119) 

The factor of sin 6* in the numerator associated with the solid angle and the sin^ in the 
denominator associated with the delta function in spherical polar coordinates cancel. The 
integrals in Eq. flll9p give 

m-^^^S{k-ko). (120) 
As for the Ising case we replace the delta function by the correlation length ^ ~ R£~^/'^. 
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16 



48 



49] and 



Hence, the structure function at k = ko scales as 

2d -1/2 

^(^o) - j^^^ = R-'e'. {d = 3) (121) 

Equation (11211) implies that there is either no divergence or the divergence is logarithmic, 
which is consistent with the fact that no direct evidence of a pseudospinodal has been 
observed in simulations of simple fluids in ci = 3, despite the indirect evidence that nucleation 
is influenced by a pseudospinodal in deeply quenched Lennard- Jones liquids 
in nickel fl7 |. 

It is easy to show that S{ko) for arbitrary dimensions scales as 

^2 D^-1/2 

~ = R-'-'e-"'-*"' « e-\ (122) 

Equation (I122p predicts that 7 = 1 for c? = 1, 7 = 1/2 for = 2, and 7 = for d = 3, in 
contrast to the mean-field result 7 = 1 for all d (see Eq. (11051) ). Results consistent with the 
predictions in Eq. (I122p were found in d = 1-3 for the potential in Eq. (I102p 50[ |. 

If we could do a measurement on a time scale of the order of the fluctuation lifetime 
e~^, we would see a smeared out density fluctuation that was radially symmetric and varied 
periodically in the radial direction. We expect that the dominant periodicity would be 
characterized by the wave vector ko and that the divergent contribution to the structure 
function near the spinodal could be approximated for k ^ k^ hj 

SUk) - I dxx'-'e^^'^-'^^. (123) 

Note that we used the density of a fluctuation rather than the cluster, which is appropriate 
for the time averaged cluster distribution. At A; = /cq the structure function will be given by 

SUk.) - ^^,R'e-''' = e-\ (124) 



consistent with the results of Refs. and 



.I45 



46 



In summary, the structure and finite lifetime of the fundamental clusters are responsible 
for the behavior of S{k) near the liquid-solid spinodal. If the peak of S{k) is at A; = as in 
the Ising model, the fact that measurements are made on a time scale short compared to the 
fundamental cluster lifetime rather than a time scale much longer as required by mean-field 
theory makes no difference to the measured value of the exponent that characterizes the 
divergence. In contrast, if the peak of S{k) is at A; 7^ 0, the measurement time scale affects 
the observed value of the exponent. 
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VIII. CLUSTER STRUCTURE AND NUCLEATION 



In classical nucleation theory [19| the nucleating droplet is assumed to be isolated, com- 
pact, and describable as a fluctuation about a quasiequilibrium metastable state. Classical 
nucleation occurs near the coexistence curve independent of the range of interaction 51|. 
For systems with sufficiently long-range interactions a quench near the pseudospinodal can 
lead to nucleation being influenced by the critical point nature of the pseudospinodal, which 
implies that the surface tension will vanish as the spinodal is approached. Near the pseudo- 
spinodal the surface tension will be nonzero, nucleation will occur with a very small surface 



tension 



ory 



10 



An e^ 
L anger 



11], and the nucleating droplet is no longer compact as in classical nucleation the- 

ll| . We refer to this form of nucleation as spinodal nucleation. 

egant way to approach nucleation theoretically was developed most fully by 



18 



19| and adapted to spinodal nucleation by Klein and Unger ll|, |52|. The 



Hamiltonian in Eq. (jl]) is used to calculate the free energy -F(/3, h) in the equilibrium state 



F{l3,h) = -fcsTln 



(125) 



which can be analytically continued from the stable to the metastable state y^. The 
nucleating droplet is associated with the solution of the Euler-Lagrange equation obtained 
from the functional derivative of the Hamiltonian |il8i] : 



i?^V>(x) - 2|e|0(x) + 40''(x) ~h = 0. 



(126) 



In Eq. (I126P e < 0, that is, the temperature T is below the critical temperature. The 
dominant exponential part of the nucleation probability is obtained by substituting the 
solution to Eq. (11261) into the expression for the imaginary part of the analytically continued 



18 and are outlined in Ref. 



19 



free energy. The details of this approach can be found in Ref. 

Klein and Unger expanded the Hamiltonian in Eq. (j4]) about the mean-field spinodal as 
in Eq. (152|) and obtained an Euler-Lagrange equation of the form 52] 



- i?W(x) + AiA/ii/2^(x) - \2^\x) = 0, 



(127) 



where Ai and A2 are constants for ixed e. It is straightforward to show that the solution to 



Eq. (11271) must have the form jll|, 



[US) 
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which imphes that the difference in the order parameter density between the interior of the 
droplet and the background is ~ A/i^/^; this difference is vanishingly small for large R and 
small Ah. This small difference presents two problems: how can we identify when and where 
the nucleating droplet occurs and how can we determine its structure. In Ising models these 
problems have been solved by mapping the spinodal onto a percolation transition and using 
the relation between the clusters and the nucleating droplet. 

Because the density of the nucleating droplet over the background is the order of A/i^/^, 
the density of the fundamental clusters is order A/i^/^/i?'^A/i^/^~'^/^ (Eq. (IHOj) ). and the den- 
sity of a fluctuation is A/i^/^/(i?'^A/i'^/^~'^/^)^/^, the nucleating droplet is not a fundamental 
cluster or a fluctuation described by the Gaussian approximation. From the discussion in 
Sec. I VI I we know that the fluctuations near a spinodal are generated by a random walk of 
the number of fundamental clusters in the up or down direction. A possible scenario might 
be that the nucleating droplet is generated by a random walk in the number of fundamental 
clusters that produces a region of density A/i^/^. From the discussion in Sees. IVl and IVll 
such a process would require a "walk" of distance /2°'A/i^/^~°'/^ {[R^Ah^^'^''^^^]'^ steps) and a 
time of r ~ A/i^^/^i?'^A/i^/^~'^/'^. Because the nucleation time is the inverse of the probabil- 
ity 13], the nucleation time would be of the order of exp{l3R'^Ah^^'^~'^^'^). Hence a random 
walk occurs too quickly to account for the time needed to see nucleation. 

A clue to the relation between the fundamental clusters and the nucleating droplets is 
provided in Ref. in which nucleation was observed near the spinodal in a = 2 Ising model 
with long-range interactions. The nucleating droplet was identified using intervention, and 
it was found that the number of fundamental clusters in a correlation length volume just 
prior to nucleation is the order of Gs, implying that Gg fundamental clusters with density 
Ah'^/'^ / gV^ coalesced into an object with a density on the order of A/i^/^. 



The results of Ref. 



121 together with the random walk argument suggests the following pic- 



ture of the relation between the fundamental clusters and the nucleating droplets. While the 
system is in the metastable state, there are fluctuations in the number of fundamental clus- 
ters. These fluctuations result in regions with order Gg = R'^Ah^^'^~'^^* fundamental clusters 
in excess of the background on a time scale ~ R'^Ah^^'^'^^^Ah^^^'^. Because this time scale is 
much less than the time scale for nucleation, which is of the order of exp(/?i?'^A/i^/^~'^/'^) 
the appearance of Gg clusters with a linear spatial extent of the correlation length will hap- 
pen everywhere in the system many times before the nucleation event. In Fig. [3] we plot 
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FIG. 3: The number of fundamental clusters in a region the size of the correlation length where 
nucleation will occur (see arrow). Note the number of times where the number of fundamental 
clusters is of order Gg- The Monte Carlo simulation was done for a d = 2 Ising model with R = 28, 
h = 1.25, Gs = 12.3, and linear dimension L = 560. 

the number of fundamental clusters the size of the correlation length in a region where we 
know that nucleation will occur. Notice the number of time intervals where the number of 
fundamental clusters is of order Gg- 

Because the time scale for nucleation is of the order of exjp((3R'^Ah^^'^~^^^) and Gs funda- 
mental clusters coalesce at nucleation, there must be a free energy cost associated with the 
coalescence that is the order of Gg. We can estimate this free energy cost by noting that 
the coalescence of Gg fundamental clusters change the density only infinitesimally. This as- 
sumption is justified because the random walk is assumed to generate a density in the stable 
phase direction of order A/i^/^ and the nucleating droplet has the same density. Therefore 
the energy change due to coalescence is negligible. The entropy cost can be estimated by 
noting that there are Gg fundamental clusters before nucleation and one nucleating droplet 
after coalescence. These considerations imply that the entropy change is given by 



~ In 2 - In 2^'^'^'''-'''' ~ -R'^Ah'/'-"/' In 2, 



(129) 



because each cluster has two states, up and down. Because the energy change is negligible, 
the free energy change due to coalescence is the order of ^ 1. Hence the probability 
of nucleation and hence the probability of coalescence is cx exp(— /3AF) ~ exp(— /^G^), in 



agreement with Refs. 



11 



and 



52 
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This discussion and the one in Sec. IVIII suggests that the fundamental clusters are not 
only a mathematical construct, but are real physical objects whose probability is not given 
by the usual Boltzmann factor (see Sec. IIIII) . This suggestion will be given further credence 
by the discussion in the next section. 



IX. CLUSTERS AND MODELS OF EARTHQUAKE FAULTS 



In this section we discuss the relevance of the cluster structure to our understanding 
of models of earthquake faults. The original Burridge-Knopoff model consists of blocks 



connected by linear springs to their nearest neighbors with spring constant [53|]. The 
blocks are also connected to a loader plate by linear springs with spring constant ki, and 
rest on a surface with a nonlinear velocity-weakening stick-slip friction force. 

A simulation is initiated by choosing the displacements of the blocks at random. While 
the loader plate is fixed, we determine the stress on each block (the force due to the springs) 
and update its velocity and displacement according to Newton's equations of motion. We 
continue these updates until all blocks are stuck. A block is "stuck" when its velocity is 



55, 



56|- We then add stress to all 



below a certain threshold and other criterion are met 
the blocks by moving the loader plate to bring the block with the largest stress to failure. 
That is, when the stress on a block exceeds the static coefficient of friction, the block "fails" 
and begins to slip. This process insures that there is only one block that initiates the failure 
sequence. An earthquake is comprised of all the blocks that fail between plate updates. 



The number Ug of earthquakes with s blocks exhibits a power law dependence on s, 



(130) 



54, 



55 



5d, 



56| . However, 



58] suggests 



with a; ~ 2 if the blocks are connected by only nearest-neighbor springs 
the calculation of realistic stress transfer Green's functions for real faults 
that we should consider springs that connect further neighbor blocks. The behavior of the 
generalized Burridge-Knopoff model is more complicated and in the limit of long-range stress 
transfer and a slow decrease of the velocity-weakening friction force with increasing velocity, 
it has been found in simulations that a; ~ 3/2 [56]. 

To provide more insight into the behavior of the model with long-range stress transfer 
we discuss a cellular automaton (CA) version of the Burridge-Knopoff model introduced 
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by Rundle, Jackson, and Brown 59|, |60|, l6l|], which considers blocks and springs as in the 
Burridge-Knopoff model. A failure threshold ap and a residual stress ctr is specified for each 
block. For simplicity, we will take the ap and an to be the same for all blocks. The stress 
on a block is given by aj = ^^(A — uj) + kc ("^j "~ "^j)) where A is the displacement of the 
loader plate and uj is the displacement of block j from its initial position. li crj < crp, we do 
nothing and proceed to the next block. If aj > ap, we move the block a distance Au where 



and q is the number of sites within the interaction range. If the range of stress transfer is 
R, then q = {2R + 1)^. Once the system is quiescent, that is, Gj < ap for all j, the plate is 
updated as in the Burridge-Knopoff model. 

The number of earthquakes with s failed blocks has a power law dependence as in 
Eq. f llSOp with x = 3/2 in the mean- field limit R ^ oo. How does this scaling arise? 
A clue is that the long-range stress transfer CA models with noise added to the stress 



63| . It also has been shown 



drop can be described by equilibrium statistical mechanics |62l . 
that this model can be described by a Langevin or Landau- Ginzburg equation in the limit 
R ^ oo [64]. The latter equation has the same form as Eq. fl52|) (see Ref.bl). These consid- 
erations imply that the scaling behavior Us ~ s~^/^ in the long-range CA models is identical 
for scaling purposes to scaling near the Ising spinodal jd]. 

To obtain the scaling exponent x, we use the fact that the Ising spinodal can be described 
by a Fisher droplet model 65| in which the system near a critical point can be described by 
a collection of non-interacting droplets. Fisher assumed that the distribution of the droplets 
scale as 16511 



-A/is* 

ns ~ (132) 

From the cluster mapping near the spinodal discussed in Appendix lAj we know that there 
exists independent objects, the fundamental clusters, which scale in a similar way. To 
calculate the exponent r for the earthquake CA model, we need only calculate the cluster 
scaling exponent for the Ising spinodal. 

To obtain the exponent r, which we will relate to the exponent x in Eq. fll30p . we note 
that the Fisher droplet model exponents are related to the spinodal exponents through the 
first several moments of Eq. f ll32p . In particular, the isothermal susceptibility x is given by 
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the second moment of hs 

-Ah. 



X oc ds oc A/^^l/^ (133) 



e 

s' 

and the order parameter density ip (see Appendix |A]) is related to the first moment of n^: 

dss^-^ (X Ah^/^. (134) 

Here we have kept Gg constant and used the fact that hyperscahng holds. This constraint 
is appropriate for tire scaling eve.ts we are corrsideriug because these models self-orga,ri.e 

to run at a fixed distance from the spinodal with Gg ~ 3-5 [3]. We can assume that the 
exponential in the integrals is approximately one until s'^ ~ Ah^^. Hence, 

3-r 1 



a 2' 

and 

r - 2 1 



(135) 
(136) 



a 2 

Equations (I135p and (11361) yield r = 5/2 which apparently differs from the measured value 
of X = 3/2. However, the exponents obtained by this reasoning assume that the cluster 
distribution is obtained by tossing bonds between occupied sites. In the earthquake case the 
clusters are grown from a seed, the site that is brought to failure by a loader plate update. 
Because for a cluster of size s there are s places that could have been the seed, the number of 

n 

such clusters is sug, where Ug is given in Eq. (11301) [39]. Hence, r = x + 1 and hence x = 3/2, 
in agreement with the simulations of models with long-range stress transfer. The scaling 
form for Ug with x ~ 2 for the usual Burridge-Knopoff model with short-range interactions 
is not understood. 

There are other properties of the distribution of earthquakes in the Rundle- Jackson- 
Brown CA model that can be obtained from consideration of the clusters in Ising models 
near spinodals. We refer the interested reader to Ref. y. The relation between the CA and 
Burridge-Knopoff models for long-range stress transfer is discussed in Refs. 



66 



and 



67 



X. SUMMARY AND CONCLUSIONS 



We have shown that the structure of the clusters near mean-field critical points and 
(pseudo) spinodals and its relation to thermal fluctuations is more complicated than the 
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corresponding relation in systems that are not mean-field and which obey hyperscaling. 
Moreover, these fundamental clusters and their structure have physical consequences which 
implies that the clusters are not only convenient mathematical constructs. 

For (near-) mean- field systems the thermal fluctuations are generated by fluctuations in 
the number of fundamental clusters. Because the clusters in (near-)mean-field systems 
play a different role and appear to be real physical objects rather than just mathematical 
constructs, we refer to them as fundamental clusters. The probability of finding a funda- 
mental cluster is not given by the Boltzmann factor because the cluster lifetime is much 
less than the decorrelation time. Their physical consequences are seen most clearly near the 
pseudospinodal in supercooled liquids where the relation between the measurement time and 
the lifetime of the fundamental clusters yields predictions for the behavior of the structure 
function that are confirmed by simulations on near-mean-field systems and are contrary to 
mean-field theory. 

In addition to the applications of the cluster structure we have discussed, there are many 
other applications that have shed light on physical processes. These applications include: 

1. The elucidation of the early time structure of systems undergoing spinodal decompo- 
sition and continuous ordering, including the understanding of why the linear theory 



of Calm, Hilliard, and Cook |20 



21 



22l | fails first at large momentum transfer 23|, 



the fractal structure of the mass distribution of early time spinodal decomposition 



24|, 



and a physical interpretation of the fermionic (Grassman) variables associated with a 



supersymmetric representation of the early stage continuous ordering 



23 



24 



68|. 



2. The phase separation of polymer and solvent in the presence of gelation 69|. 

3. Possible precursors to nucleation near the pseudospinodal [l^ . 

Future work includes the application of cluster methods to the study of precursors to large 
earthquakes in the CA and Burridge-Knopoff models, the study of heterogeneous nucleation 
near pseudospinodals, the investigation of fracture and the merging of microcracks, and 
the study of the crossover from the linear regime of spinodal decomposition and continuous 
ordering to nonlinear evolution. 
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APPENDIX A: PERCOLATION MAPPING 

To obtain a deeper understanding of the structure of the fluctuations near the mean- 
field critical point, we map the Ising critical point (mean-field and non-mean-field) onto a 



percolation transition for a properly chosen percolation model 



We first describe the 



mapping introduced by Kasteleyn and Fortuin [70|] of the s-state Potts model onto random 
bond percolation. The latter is defined on a lattice where all the sites or vertices are occupied 
with probability one, and the bonds are occupied with a probability pf,. Clusters are defined 
as a set of sites connected to each other by bonds and not connected to any other sites in 



the lattice [39 1 



The Hamiltonian for the s-state Potts model is 

Hp = -Jp 5^(5.,., - 1) - /^P(5.a - 1), (Al) 

where (Tj specifies the state of site i, Jp > is the coupling constant, and hp is the Potts 
field; the Kronecker delta S^j.^^ ^ only when sites i and j within the interaction range are 
in the same state. We first set hp = and write the Boltzmann factor e~^^^ as 

e-^^^ = n K.^. + - ^...)] , = n [(1 - --''nS..^. + ^-''1 • (A2) 

ij ij 

We associate a bond with Sa^aj = 1 and the absence of a bond with 5^^^^ = 0. With this 
association the generating function for the random bond percolation model is obtained by 
differentiatirig the free energy for the s-state Potts model with respect to s and then setting 
s equal to 1 [70] . There are many subtle mathematical points in the Kasteleyn- Fortuin proof 



70 



and the references therein for the details. 



of this relation, and we refer the reader to Ref. 
Because it will be needed to understand the structure of fluctuations, we demonstrate how 
the connection between the Potts model and percolation works. 
The partition function Zp is given by 

Zp = ^e~^^p. (A3) 

and the free energy in the canonical ensemble is Fp{f3, s) = —ksT In Zp. If we differentiate 
— /3Fp(/5, s) with respect to s, we obtain 

P ds -l^Zp ds ■ ^^^^ 
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Setting s = 1 results in Hp = because there is only one Potts state; hence Zp = = 1 
for s = 1, where is the number of sites in the lattice. Therefore 

_ dFp^ =^ye-^-^ . (A5) 

OS s=l OS ^ s=\ 
c 

The percolation generating function is the right-hand side of Eq. (lASp . To understand 
this interpretation we consider several terms in Eq. (jASj. We use Eq. (IA2I) for e~^^^ and 
first consider the term e'^"^^ in each of the factors in the product; that is, we include no 
terms with . For a lattice with c = qN/2 total possible bonds, we find a contribution 
to Fp{f3,s) of the form s^e"^'^^'^^. By differentiating with respect to s and setting s = 1, 
we obtain the contribution to the generating function Gf. 

Gf,i = Ne-f^-^^"^. (A6) 

Because there are A^ sites and cN bonds, Gf^i can be interpreted as the mean number of 
single site clusters. That is, e~^''^'^^ is the probability that there are no bonds present. 

We now consider a term from Eq. (lA2p that includes only one delta function, which we 
take to be 5(jio-2- The contribution to Gf has the form 

Gf,2(p) = (1 - e-^^p)se-^-^p(^^-i)s^-2 = ^^(1 - PkT''-h''-\ (A7) 

where we have associated the bond probability with l—e~^^^. Differentiating with respect 
to s and setting s = 1 gives A^ — 1 for the number of clusters times the probability of such 
a configuration. There are N — 2 one-site clusters and one two-site cluster for a particular 
bond. The number of ways we can choose one bond is cN, so that the first two contributions 
to dF{(3, s)/ds are 

Gi,iiP) + GiM = Nil-pbT'' + {N- l)cNp,{l-p,r''-\ (A8) 

If we continue in this manner, we would find that the terms we obtain are the number of 
clusters in a given configuration. The complete enumeration of the configurations will lead 
to an expression for the mean number of clusters as a function of pb. Hence, the mean 
number of clusters can be written as 

G'f(p) = 5^<nfc>, (A9) 
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where <nk> is the mean number of clusters with k sites. 

To obtain the full generating function for random bond percolation, we must include the 
field hp in the calculation of the free energy. We write 

,-0H, = [(1 _ e-^'nS..., + e-^'-] n [(1 - e-f^'^nS.a + e"^'^^] . (AlO) 

ij I 

The terms in the expansion of Eq. (lAlOp represent sets of connected sites (clusters) generated 
by Terms of the form S^j^i, where the index 1 labels one of the s possible states of a 

site, give no contribution to the derivative of the partition function with respect to s because 
the 6cr^.i term fixes all spins in a cluster to the Potts state labeled as 1. Hence there is no 
s dependence and no factor of s in the product, which equals s raised to the power of the 
number of clusters. From Eq. ( lAlOp all clusters with a nonzero weight after differentiation 
with respect to s will have a field dependence of the form e~^^^, where n is the number of 
sites in the cluster. If we resum as in Eq. ( lASp . we obtain [3] 



Gf (Pfc, hp) = J2 <nk> e-^"''^. (All) 

k 

The reader might want to work out the generating function for small lattices to see how the 
sum in Eq. (I All!) arises. 

To investigate the mapping of the percolation model onto the Ising model, we consider 
the dilute s-state Potts [71] model with the Hamiltonian 

PH^P = -(3Jp - iKn, - (3hp 5^(5,a - 1)^* - ^LG n^rij + A ^ n^. (A12) 

ij i ij i 

where Ui = 1 denotes that a particle (spin) occupies site i and nj = denotes an empty site. 
Hence there is a Potts interaction between occupied (n^ = 1) sites. The quantity 

ij i 

is the Hamiltonian for the lattice gas formulation of the Ising model \(\, with Ki^q the 
(dimensionless) coupling constant and A the chemical potential. In terms of the parameters 
in the Ising Hamiltonian Hi, 

ij i 
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with Sj = ±1, we have 72 1 



i^LG = 47^1 and A = phi + 2cKi. 



(A15) 



Differentiating the free energy constructed from iJ^p with respect to s and setting s = 1 
results in the generating function for correlated site random bond percolation for which oc- 
cupied sites are distributed according to the lattice gas Hamiltonian in Eq. (]A13I) and bonds 
are thrown randomly with a probability pb between pairs of occupied sites. To understand 
this connection we note that 



^DP 



(A16) 



and 



ds 



ds (2^{ffi}{ni} ^ 



s=l 



(A17) 



^DP 



If we write 



and 



exp [(3Jp{6^^^^ - l)n,nj] = [(1 - e-^^^)6^^^^ + e-^-^^]ninj + (1 - riirij), 



exp [Php J2{6^a - = (1 - e-^'^n^.^mi + e-^'^n, + (1 - n,), 



(Al8a) 



(Al8b) 



we can use the same arguments that we gave for random percolation to show that the 
expression in Eq. (1A17I1 leads to the generating function for correlated site random bond 
percolation where the occupied sites are distributed according to the lattice gas Boltzmann 
factor constructed from the Hamiltonian in Eq. flAlSP 711 ]. Note that for this model the 
sum over the s states of the Potts spin is one for all s if site i is empty. 

We now consider the Hamiltonian H^p in Eq. flA12p . We set hp = hj = 0, and hence 
A = 2qKi = qKi^c/^ from Eq. fIXTSll and let Jp = in Eq. flXT2l) . Then Hp^p becomes 



DP 



(A19) 



ij ij ij 

Suppose that for a pair of sites within the interaction range, either both sites are empty or 
both sites are filled, but the spins are in the same Potts state. When both sites are empty, 
there is only one Potts configuration. When both sites are occupied and in the same Potts 
state, there are s configurations. Hence, there are s + 1 configurations with Hp,p = 0. Now 
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consider the case where either one site is occupied and one empty or both sites are occupied, 
but the spins are in different Potts states. In this case the contribution of this pair to H^p 
is Ki^Q,/2, and the number of ways this combination can be obtained is (s + l)s. These 
considerations imply that the dilute s-state Potts model at Jp = Klg/2 is equivalent to a 
pure (s + l)-state Potts model with the Hamiltonian 

/?i/p,(.+i) = -^(5.,.,-l), (A20) 



where can be in s + 1 states. The s = 2 Potts model is the lattice gas model in Eq. ( 1A12I1 



with A = cKi^q/2. That is, for (3Jp = Ki^q/2 = 2Ki and hj = hp = 0, the Hamiltonian of 
the dilute s-state Potts model is the same as the Hamiltonian of the {s + l)-state pure Potts 
model. In the limit s — 1, the (s + l)-state pure Potts model is the lattice gas model. 
In this formulation we can write the Boltzmann factor as 



g-/3/^P,(.+i) = - ^ (A21) 

Clearly the singular behavior of the free energy comes from the terms in Eq. flA22p contained 
in the product over lattice sites. This product has the form of Ylij [fij + 1] ? where we associate 
a graph or cluster with a product of the fij summed over Potts states. These clusters are 
the same as the percolation clusters because the sites are connected by 5o- a- bonds. The 

nn 

linked cluster theorem [73|, [7^ states that the singular part of the free energy Fp^sing is the 
sum over all connected graphs in the thermodynamic limit. Connected graphs are those in 
which all points or vertices of the graph are connected by an fij bond. Because the graphs 
are connected, the sum over Potts states results in a factor of s, independent of the size or 
structure of the graph. Therefore the derivative of -Fp sing with respect to s results in the 
same sum without the overall factor of s. Thus 

which implies that the percolation transition and the Ising critical point occur at the same 
temperature and have the same critical exponents. The amplitudes of the singular quantities 
differ by a factor of two. If instead of the Hamiltonian in Eq. (]A12p . which defines clusters 
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as consisting of only occupied sites, we add a term to the Hamiltonian of the form 

/3i?DP, empty = -{^Jv " " " ^j) 

- i^LG 5^(1 - (1 - ^j) + A' 5^(1 - nO, (A24) 

ij i 

we can define clusters between "empty" sites. In this way the singular part of the free energy 
will be identical to the mean number of clusters for s — > 1 and /ip = 0. The same result was 



obtained using renormalization group techniques in Ref. 138|. The mapping of the Ising model 
onto a properly chosen percolation model is not restricted to nearest neighbor interactions. 
If the Ising interaction has a range i?, we need only to choose the Potts interaction to also 
have a range i?, which means that the bond between sites is randomly placed between any 
two spins within the interaction range. 

We now consider the mapping of a thermal problem onto a percolation model near the 
spinodal. This mapping will require a slightly different approach. We again begin with the 
dilute s-state Potts model. The Hamiltonian is the sum of ifop and i^LG from Eqs. ( ]A12I) and 
( IA13I1 . This Hamiltonian can be put into a continuum form using the Gaussian transforma- 
311 . 175l | . Because we are interested in the mean-field limit , the Landau-Ginzburg- Wilson 

and 



25 



75 



tion 

Hamiltonian is the free energy. We have from Refs 

i^Dp(C, 0) =Jd^[^s{s - l)(i?VC(x))' - ns{s - l)C'(x) - hp{s - l)C(x) 

+ ^s{s -l)is- 2)is - 3)C^(x) + ^s{s - l)C'(x)0(x)] + (A25) 

where is the free energy in Eq. (HHj) . The constants ri, wi, W2, and e in Eq. PHI) can 
be written as functions of J, Ki^q, and c. The global percolation order parameter ( is the 
probability that a spin in the stable phase direction belongs to the infinite cluster of occupied 
sites. As for a discrete system, the percolation model is obtained by differentiating F^p with 
respect to s and setting s = 1 3]. 

To map the thermal problem near the spinodal onto the percolation problem, we rewrite 
Eq. pHl) as in Eq. fl52l) with Ah = hg — hj. We equate the functional derivative of 
dFj)p{(, (f))/ds\s=i = Fp with respect to C(x) to the functional derivative of F{ip) with 
respect to ip- That is, 

^ = -i?V^C(x) - 2riC(x) + ^C'(x) + ^2C(x)</.(x) - hp, (A26) 

and 
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= -i?2vV(x) + 2AiA/ii/V(x) - 3A2^'(x) + 4A3^'(x) (A27) 

dip 

must be equal. This condition implies that — 2(ri — ^0s)C + ^ must be the same 
as 2AiA/i^/^^/' — 3A2'?/'^. We dropped the term in the Euler-Lagrange equation because 
as we shown in Sees. VI and XI. 
At the spinodal {Ah = 0) we identify ( with —ip and require that hp = 0, Wi/3 = — 3A2, 
and ri = W20(x)/2. With these equahties the solutions of Eqs. (1A26P and (1A27|) with 



5Fp/5^(x) = 6F{ip)/6ip{'x) = are identical. If we write the parameters in Eqs. flA26P and 
(IA27P in terms of the parameters Jp, Ki^g, and c of i^DP in Eqs. (IA12p and (lAlSp . we obtain 

pjp = 2Kj{l-p), (A28) 

where the density p = (1 + m)/2 and (p is proportional to m, the magnetization per spin. 
Because we are interested in the coincidence of the spinodal with a percolation transition, 
we can set 0(x) = 0^, the value of the order parameter at the spinodal. 

To define the bond probability in terms of a physically meaningful quantity, we need to 
obtain the proportionality factor between (ps and m^, the value of the magnetization at the 
spinodal. To do so we use the relation between the Ising model and the 0^ theory generated 
by the Hubbard- Stratonovich transformation [31|. In particular, the value of (ps is [3] 

When T or equivalently Kj —>■ 00, rris ^ 1 as can be seen by noting that the magnetic 

1 /2 

field is divided by T in the Boltzmann factor. Because (ps — > 1/cKj as Kj 00, we have 

m, = cKy^(ps, (A30) 



so that m-s ^ 1 as Kj 00. If we use Eq. flASOp . we obtain the expression for the bond 



probability that maps the percolation model onto the spinodal 

Pb=l-e-^■^^ (A31) 

where Jp is given in Eq. flA28p . The validity of this mapping was demonstrated numerically 
in Refs.HQ, andQ 



The interpretation of this mapping is that the spinodal curve is the locus of a set of 
percolation transitions. If the spinodal curve is approached with h 0, there is a transition 
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to a spanning cluster in the stable phase direction at the spinodal. We stress that this result 
is correct only in the mean-field limit, — > oo. 

The various mappings we have discussed imply that the free energy of the lattice gas 
model is isomorphic to the generating function for correlated site random bond percolation. 
For hp = hi = 0, the generating function is the mean number of clusters. 
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